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Abstract 

Finite N effects on the time evolution of fuzzy 2-spheres moving in flat spacetime are studied 
using the non-Abelian DBI action for N DO-branes. Constancy of the speed of light leads 
to a definition of the physical radius in terms of symmetrised traces of large powers of 
Lie algebra generators. These traces, which determine the dynamics at finite N, have a 
surprisingly simple form. The energy function is given by a quotient of a free multi-particle 
system, where the dynamics of the individual particles are related by a simple scaling of 
space and time. We show that exotic bounces of the kind seen in the 1/N expansion do 
not exist at finite N. The dependence of the time of collapse on N is not monotonic. The 
time-dependent brane acts as a source for gravity which, in a region of parameter space, 
violates the dominant energy condition. We find regimes, involving both slowly collapsing 
and rapidly collapsing branes, where higher derivative corrections to the DBI action can be 
neglected. We propose some generalised symmetrised trace formulae for higher dimensional 
fuzzy spheres and observe an application to D-brane charge calculations. 
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1 Introduction 


The symmetrised trace prescription for the non-Abelian action of multiple DO-branes was 
proposed in j[Q and extended to include background RR fluxes in [2]. An interesting time 
dependent system, in which the need for an exact prescription arises, is a spherical bound 
state of N DO-branes with a spherical D2-brane, for finite values of N. This can be studied 
both from the point of view of the Abelian D 2 DBI action and the non-Abelian DO-DBI 
action. The latter configuration also has an M-theory analogue, that of a time dependent 
spherical M2-brane, which has been studied in the context of matrix theory EH 03- In |S] 
it was shown that the DO-brane construction, based on the fuzzy 2-sphere, agrees with the 
Abelian D2-construction at large N. 1/N corrections coming from the symmetrised trace 
and a finite N example were also studied. Here we develop further the study of finite N. 
The need for the non-linear DBI action as opposed to the Yang-Mills limit of the lower 
dimensional brane was recognised in a spatial D1 _L D 3 analog of the DO — D 2 system ||. 

In this paper we extend the calculation of symmetrised traces from the spin half example 
of 0 to general representations of SO(3). These results allow us to study in detail the 
finite N physics of the time-dependent fuzzy two-sphere. We begin our finite N analysis 
with a careful discussion on how to extract the physical radius from the matrices of the 
non-Abelian ansatz. The standard formula used in the Myers effect is R 2 = Tr($i$j)/IV. 
Requiring consistency with a constant speed of light, independent of N, leads us to propose an 
equation in section 2, which agrees with the standard formula in large N commutative limits, 
but disagrees in general. Section 3 gives finite N formulae for the energy and Lagrangian of 
the time-dependent fuzzy 2-sphere. We also give the conserved pressure which is relevant for 
the Dl _L D 3 system. In section 4, we study the time of collapse as a function of N. In the 
region of large IV, for fixed initial radius Rq, the time decreases as N decreases. However, 
at some point there is a turn-around in this trend and the time of collapse for spin half is 
actually larger than at large N. We also investigate the quantity E 2 — p 2 , where E is the 
energy and p the momentum. This quantity is of interest when we view the time-dependent 
D-brane as a source for spacetime fields. E is the T 00 component of the stress tensor, and 
p is the T 0r component as we show by a generalisation of arguments previously used in the 
context of BFSS matrix theory. For the large N formulae, E 2 —p 2 is always positive. At finite 
N, this can be negative, although the speed of radial motion is less than the speed of light. 
Given the relation to the stress tensor, we can interpret this as a violation of the dominant 
energy condition. The other object of interest is the proper acceleration along the trajectory 
of a collapsing D2-brane. We find analytic and numerical evidence that there are regions of 
both large R and small R , with small and relativistic velocities respectively, where the proper 
accelerations can be small. This is intriguing since the introduction of stringy and higher 


2 


derivative effects in the small velocity region can be done with an adiabatic approximation, 
but it is interesting to consider approximation methods for the relativistic region. 


In section 5, we discuss the higher fuzzy sphere case 0 0 ei uni im ca CHI d • We give a 
general formula for STr(XiXi) m , in general irreducible, representations of SO(2k + 1). This 
formula is motivated by some considerations surrounding D-brane charges and the ADHM 
construction, which are discussed in more detail in |30j . Some of the motivation is explained 
in Appendix A. This allows us a partial discussion of finite N effects for higher fuzzy spheres. 
We are able to calculate the physical radius following the argument of section 2; however, in 
general one needs other symmetrised traces involving elements of the Lie algebra so(2k + 1). 

The symmetrised trace prescription, which we study in detail in this paper, is known to 
correctly match open string calculations up to the first two orders in an a! expansion, but the 
correct answer deviates from the (c/) 3 term onwards deni Emm- It is possible however 
that for certain special symmetric background configurations, it may give the correct physics 
to all orders. The D-brane charge computation discussed in the Appendix can be viewed as 
a possible indication in this direction. In any case, it is important to study the corrections 
coming from this prescription to all orders, in order to be able to systematically modify it, if 
that becomes necessary when the correct non-Abclian D-brane action is known. Conversely 
the physics of collapsing D-branes can be used to constrain the form of the non-Abclian 


DBI. 


2 Lorentz invariance and the physical radius 


We will study the collapse of a cluster of N IAO-branes in the shape of a fuzzy S 2k , in a 


flat background. This configuration is known to have a large-A" dual description in terms 
of spherical D(2k) branes with N units of flux. The microscopic DO description can be 


obtained from the non-Abelian action for a number of coincident branes, proposed in 



( 2 . 1 ) 


where 



( 2 . 2 ) 


Here a, b are worldvolume indices, the c&’s are worldvolume scalars, A = 2ir£ 2 and 


Qij — Sij + . 

We will consider the time dependent ansatz 


(2.3) 


= R(t)Xi , 


(2.4) 
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The Xi are matrices obeying some algebra. The part of the action that depends purely on 
the time derivatives and survives when R = 0 is 


Sdo — 


f dtSTry/1 - \ 2 {d t ^y 


dtSTr\J 1 - A 2 (d t R) 2 X i X i . 


(2.5) 


For the fuzzy S 2 , the X ^ = cq, for i = 1,2,3, are generators of the irreducible spin n/2 
matrix representation of su( 2), with matrices of size N = n + 1. In this case the algebra is 


[cr^, o^] 2ieijkOLk 

and following [5], the action for N HO-branes can be reduced to 

s„ = --h 

9s^s 


dt STry 1 + 4A 2 R‘ i aia i y 1 — \ 2 (d t R) 2 onoii . 


If we define the physical radius using 


Ri 


= V Hm flhhhry = lim 


( 2 . 6 ) 


(2.7) 


( 2 . 8 ) 


phys m—>oo m—»oo ST r (cK,Q(j) r 

we will find that the Lagrangian will be convergent for speeds between 0 and 1. The radius 
of convergence will be exactly one - this follows by applying the ratio test to the series 
expansion of 

STr\Jl- X&onon , (2.9) 

where a dot indicates differentiation with respect to time. This leads to 


R%„ = A 2 RV 


( 2 . 10 ) 


Using explicit formulae for the symmetrised traces we will also see that, with this definition 
of the physical radius, the formulae for the Lagrangian and energy will have a first singularity 
at Rph.ys = 1- In the large n limit, the definition of physical radius in (12.101) agrees with 
.2], where R p h ys is defined by R 2 hys = Note that this definition of the physical 

radius will also be valid for the higher dimensional fuzzy spheres, and more generally in any 
matrix construction, where the terms in the non-Abelian DBI action depending purely on 
the velocity, are of the form \J 1 — X 2 X i X i (d t R) 2 . 

In what follows, the sums we get in expanding the square root are conveniently written 
in terms of r, s, defined by r 4 = 4A 2 R 4 and s 2 = A 2 R 2 . It is also useful to define 

= 

2 ’ 


f 2 = 


d2 

^phys 9 

£ 2 = r 71, 


*2 2 2 
s — s n 


( 2 . 11 ) 


The r and s variables approach the variables called r, s in the large n discussion of [^j. 
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3 The fuzzy S 2 at finite n 


For the fuzzy S 2 , the relevant algebra is that of su( 2), equation (|2.(ij) above. We also have 
the Casimir 

c = cqcq = ( N 2 — 1) , 

where the last expression gives the value of the Casimir in the TV-dimensional representation 
where N — n + 1, and n is related to the spin J by n = 2 J. 

We present here the result of the full evaluation of the symmetrised trace for odd n 

1 x (n+l)/2 

J E (a-ir, (3.i) 

i=l 

whilst for even n 

Th 

C(m, n) = -^—STria^r = 2{ ? m + 1} V(2z) m . (3.2) 

n + 1 7i+l ' 

i =1 

For m = 0 the second expression doesn’t have a correct analytic continuation and we will 
impose the value STr(aiai) () = 1. The expression for C(m, 1) was proved in J2|. A proof of 
(HOI) for 77 = 2 is given in Appendix B. The general formulae given above are conjectured on 
the basis of various examples, together with arguments related to Z+brane charges. These 
are given in Appendix A. There is also a generalisation to the case of higher dimensional 
fuzzy spheres, described in section 6 and the Appendices. 

We will now use the results (ED), (E3), to obtain the symmetrised trace corrected energy 
for a configuration of N time dependent ZlO-branes blown up to a fuzzy S 2 . The reduced 
action (EH) can be expanded to give 


C(m,n ) =- STr(aiai) m = 

77+1 


2(2777 + 
77+1 


£ 


-STr 1 + 4A 2 A 4 cqcq \/1 - A 2 k 2 a i a i 
—STr \/1 + r 4 ai.oti \/l — s 2 a.iOii 


-STr ,r + ' 

m =0 Z=0 




(- 1 ) 


( 3 . 3 ) 

( 3 . 4 ) 


The expression for the energy then follows directly - 

£ = -STr r Al (2m - 1 )(a i a i ) m+l f 1 ^ (- 1 )”\ ( 3 - 5 ) 

m =0 1=0 V / V / 

and after applying the symmetrised trace results given above we get the hnite -77 corrected 
energy for any finite-dimensional irreducible representation of spin-^ for the fuzzy S 2 . 
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For n — 1,2 one finds 


1 . . 1 + 2 r 4 — r 4 s 2 

-£ n =i(r,s) = 

l c ^ _ 2 (1 + 8r 4 - 16r 4 s 2 ) | 1 

3 n=2(r ’ “ 3 a/1 + 4r 4 (l — 4s 2 ) 3 / 2 + 3 ' 


(3.6) 

(3.7) 


We note that both of these expressions provide equations of motion which are solvable by 
solutions of the form r = t. 


For the case of general n, it can be checked that the energy can be written 

2 - 2(2 1 - l)V((2/ - 1)V - 2) 


n+1 

2 


£n(r,s) = 


y/1 + (21 - l) 2 r 4 (l - (21 - 1 ) 2 S 2) 3 /2 ’ 


for n-odd, while for n even 


S'r 4-U2T 1 - 8fr4 < 2 ^- 1 ) 

n{ ’ 1 Vl + 4(V(1 - 4 Ps^ ' 

Equivalently, the closed form expression for the Lagrangian for n odd is 

1 — 2 l 2 s 2 + / 2 r 4 ( 2 — 3/ 2 s 2 ) 


n+1 

2 


£ n (r,s) = — 2 , _ __ 

a/1 + / 2 r 4 v / l — l 2 s 2 


whilst for 77 -even 


£ n (r,s) = -1 - 2^ 


1 - 2/ 2 s 2 + l 2 r 4 (2 - 3l 2 s 2 ) 


\/l + / J r' v t — 


(3.8) 


(3.9) 


(3.10) 


(3.11) 


It is clear from these expressions that the equations of motion in the higher spin case will 
also admit the r — t solution. Note that, after the rescaling to physical variables dim 
these Lagrangians have no singularity for fixed r, in the region 0 < s < 1. In this sense they 
are consistent with a fixed speed of light. However, they do not involve, for fixed r, the form 
a Jdt 2 — dr 2 and hence do not have an so(l, 1) symmetry. It will be interesting to see if there 
are generalisations of so(l, 1), possibly involving non-linear transformations of dt,dr , which 
can be viewed as symmetries. 

One can also get exact results for the symmetrised trace corrected pressure of the fuzzy- 
S 2 funnel configuration. The relationship between the time dependent HO — D 2 system and 
the static Dl _!_ D 3 intersection was established in ESI- In that paper, the large-77 behaviour 
of both systems was described by a genus one Riemann surface, which is a fixed orbit in 
complexified phase space. This was done by considering the conserved energy and pressure 
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and complexifying the variables r and dr = s respectively. Conservation of the energy- 
momentum tensor then yielded elliptic curves in r, s, involving a fixed parameter r 0 , which 
corresponded to the initial radius of the configuration. For our system we simply display the 
general result and the first two explicit cases 


V 


^V n= i(r,s) 
\v n =2{r, s ) 


STr ± ± S 2 ”V- !)(«,«,)”+' (4 (f 

m =0 1=0 ' ' ' 

1 + 2 r 4 + r 4 s 2 
\/l + r 4 (l + s 2 ) 3 / 2 
2(1 + 8r 4 + 16r 4 s 2 ) 1 

_ 3 +1 + 4r 4 (l + 4s 2 ) 3 / 2 ~~ 3 ' 


(3.12) 

(3.13) 

(3.14) 


Similar results to those for the time dependent case hold for the exact expression of the 
pressure for the general spin-1) representation. Note again that these expressions will provide 
equations of motion which are solved by solutions of the form r = 1/cr, where a is the spatial 
D 1 worldvolume coordinate. An easy way to see this is to substitute s 2 = r 4 in turn ra, 
to find that the pressures become independent of r and s. Since the higher spin results for 
the pressure are sums of the n — 1 or n — 2 cases, the argument extends. 


3.1 Finite N dynamics as a quotient of free multi-particle dynamics 

Using the formulae above, we can see that the fuzzy S 2 energy for general n is determined 
by the energy at n — 1. In the odd n case 

n+1 n+1 

C(m, n) = +_ C(m, 1) T(2t 3 - l) 2 ” = —'-(2m + 1 ) +(2 !3 - l) 2 ” . 
n+1 n + 1 

*3 = 1 *3 — 1 

Using this form for C(m,n ) in the derivation of the energy, we get 

n+1 

2 

£Jr, s) = ^ £ n=1 (ry/(2i 3 - 1 ) , s(2i 3 - 1 )) . ( 3 . 15 ) 

*3 = 1 

Similarly, in the even n case, we hnd 

n 

2 

£n(r,s) = y £ n =‘i{;r\pr 3 , s(i 3 )) . ( 3 . 16 ) 

*3 = 1 


It is also possible to write C(m, 2) in terms of C(m, 1) as (for m ^ 0) 


22m+l 22^+1 

C(m, 2) = ———C(m, 1) = —(2m + 1) . 


(3.17) 
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Thus we can write £ n (r, s ), for even n, in terms of the basic £ n= \ (r, s) as 


£ n (r,s) = 1 + 5^ £n=i{r\/{2i 3 ) , s(2i 3 )) . (3.18) 

*3 = 1 

These expressions for the energy of spin n/2 can be viewed as giving the energy in terms 
of a quotient of a multi-particle system, where the individual particles are associated with 
the spin half system. For example, the energy function for (n + l)/2 free particles with 
dispersion relation determined by £ n =i is JT £ n= 1 (rj, s*). By constraining the particles by 
r t = r\j2i + 1 , s* = s{2i + 1 ) we recover precisely (ESP - 

We can now use this result to resolve a question raised by |j> on the exotic bounces 
seen in the Lagrangians obtained by keeping a finite number of terms in the 1 jn expansion. 
With the first 1/n correction kept, the bounce appeared for a class of paths involving high 
velocities with 7 = - 7=5 ~ c 1 / 4 , near the limit of validity of the 1/n expansion. The bounce 
disappeared when two orders in the expansion were kept. It was clear that whether the 
bounces actually happened or not could only be determined by finite n calculations. These 
exotic bounces would be apparent in constant energy contour plots for r, s as a zero in the 
first derivative dr/ds. In terms of the energies, this translates into the presence of a zero 
of dE/ds for constant r. It is easy to show from the explicit forms of the energies that 
these quantities are strictly positive for n — 1 and n — 2 . Since the energy for every n can 
be written in terms of these, we conclude that there are no bounces for any finite n. This 
resolves the question raised in |>j about the fate at finite n of these bounces. 

We note that the large-n limit of the formula for the energy provides us with a consistency 
check. In the large n-limit the sums above become integrals. For the odd-n case (even-n can 
be treated in a similar fashion), define x = 28,3-1 ~ —. Then the sum in (13.1511 goes over to 
the integral 

n f 1 2 — 2 x 2 n 2 r 4 (x 2 n 2 s 2 — 2 ) nV 1 + r 4 n 2 . 

— / ax , -- = — , . 3.19 

2 J 0 \/l + x 2 n 2 r A { 1 — x 2 s 2 ) 3 / 2 \/l — s 2 n 2 

By switching to the f, s parameters the energy can be written as /Y 1 /"/) 4 . This matches 
exactly the large n limit used in [5j. 

4 Physical properties of the finite N solutions 

4.1 Special limits where finite n and large n formulae agree 

In the above we compared the finite n formula with the large n limit. Here we consider 
the comparison between the fixed n formula and the large n one in some other limits. On 

















physical grounds we expect some agreement. The DO — D 2 system at large r and small 
velocity s is expected to be correctly described by the D 2 equations. These coincide with 
the large n limit of the .DO. In the D 1 _L D 3 system, the large r limit with large imaginary 
s is also described by the D 3. 

Such an argument should extend to the finite-n case. In ra. these systems were simply 
described by a genus one Riemann surface. However, in this case the energy functions are 
more complicated and the resulting Riemann surfaces are of higher genus. We still expect 
the region of the finite n curve, with large r and small, real s, to agree with the same limit 
of the large n curve. We also expect the region of large r and large imaginary s to agree 
with large n. 

For concreteness consider odd n. Indeed for large r, small s, (El gives 

^^4(2/ — l)r 2 ~ nr 2 , (4.1) 

i 

which agrees with r 2 . In this limit, both the genus one curve and the high genus finite n 
curves degenerate to a pair of points. Now consider large r and large imaginary s. This is 
the right regime since the D 1 _L D 3 system is described by r ~ - which means that r is 
large at small a, where ^ = is is large. For s = iS 

V ~ -nr 2 IS , (4.2) 

which agrees with the same limit of the large n curve. In this limit, both the large n genus 
one curve and the finite n curves of large genus degenerate to a genus zero curve. 

The agreement in El between the DO and D 2 pictures is a stringy phenomenon. It 
follows from the fact that there is really one system, a bound state of DO and D 2 branes. A 
boundary conformal held theory would have boundary conditions that encode the presence 
of both the HO and D2. In the large N limit, the equations of motion coming from the 
.DO-effective action agree with the H2-effective action description at all R. This is because 
at large N it is possible to specify a DBI-scaling where the regime of validity of both the HO 
and H2 effective actions extends for all R. This follows because the DBI scaling has £ s —> 0 
|20j . Indeed it is easy to see that the effective open string metric discussed in EDI has the 
property that £ 2 S G~ X = i goes to zero when N —> oo with L = £ s Vtt iV, R fixed. This 
factor £ 2 G -1 controls higher derivative corrections for the open string degrees of freedom. 
At finite N, we can keep £ 2 s G~ l small, either when R <C L or R L. Therefore, there are 
two regimes where the stringy description reduces to an effective field theory, where higher 
derivatives can be neglected. The agreement holds for specified regions of R as well as s, 
because the requirement £ 2 G -1 -C 1 is not the only condition needed to ensure that higher 
derivatives can be neglected. We also require that the proper acceleration is small. At large 
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R , the magnetic flux density is small (as well as the higher derivatives being small) and 
the 142-brane without non-commutativity is a good description. This is why the finite N 
equations derived from the HO-brane effective field theory agree with the Abelian D2-picture. 
For small R , small s, we can also neglect higher derivatives. This is the region where the 
DO- Yang-Mills description is valid, or equivalently a strongly non-commutative _D2-picture. 


4.2 Finite N effects : Time of collapse, proper accelerations and violation of 
the dominant energy condition 

We will consider the time of collapse as a function of n using the definition of the physical 
radius given in section (2). In order to facilitate comparison with the large n system, we will 
be using f, s variables. To begin with, consider the dimensionless acceleration, which can be 
expressed as 

f),f. I, 

(4.3) 


—s 


df£ I 


with 7 = 1 / yj(l — s 2 ). As the sphere starts collapsing from r = f 0 down to f = 0 , the speed 
changes from s = 0 to a value less than s — 1. It is easy to see that the acceleration does 
not change sign in this region. Using the basic energy S = S/N from fKl.(il) . we can write 


d£n=i (r,S) 
ds 

d£ n =i(r,s) 

dr 


= s- 


. (3(1 + f 4 ) + f 4 (l — s 2 )) 
V(l + r 4 )(l-s 2 )l : 

2 f 3 


(1 + r) + (1 — s J )(2 + r 4 ) 


(4.4) 


(1 + r 4 ) 2 (1 — s 2 ) 2 

Neither of the partial derivatives change sign in the range s — 0 to 1. Hence the speed s 
increases monotonically. The same result is true for n > 1, since the energy functions for all 
these cases can be written as a sum of the energies at n = 1 . 

In the n — 1 case ,r — r, s — s. For fixed ro the speed at r = 0 is given by 


(1 


„ 2 | x _ ( 1 + r o ) 5 

* 71=1 ) , s 2 ' 

(l + 2r‘)i 

Comparing this with the large n formula 

(1 - S 2 |n=oo) = (1 +ro) _1 , 

it is easy to see that 


(1 - S 2 )|„, 0 


> 2 )U ) (l + r„ 4 ) 4<1 ’ 

which establishes that the speed at r = 0 is larger for n — oo. 


(1 + 2 r ‘) 2 


(4.5) 


(4.6) 


(4.7) 
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We can strengthen this result to show that the speed of collapse at all r < ro is smaller 
for n — 1 than at n = oo. For any r < r 0 we evaluate this energy function with the speed of 

4_ 4 

collapse evaluated at s 2 = r \ , r 1 , which is the speed at the same r in the large n problem. Let 

UH 1 _ 


us define F(r, r 0 ) = £ n =\ [ r, s — 


We compare this with £ n =i(r, s ) for s appropriate 


for the n — 1 problem, which is just 1 Ldin = G(ro) by conservation of energy. We now use 

v 1+r o 

the fact, established above, that ^‘y 1 is positive for any real r. This means that we can 


show s\ n= i < 



by showing that F(r, r 0 ) > G(r 0 ). A short calculation gives 


F(r,r 0 ) -G(r 0 ) 


v 1 + 'u!'l + r 4 ) 



(4.8) 


It is clear that we have the desired inequality, showing that, at each r, the speed s in the 
n — 1 problem is smaller than the speed in the n = oo system. Hence the time of collapse 
is larger at n — 1. In the n = 2 case, we find that an exactly equivalent treatment proves 
again that the collapse is slower than at large n. However, this trend is not a general feature 
for all n. In the leading large -N limit, the time of collapse is given by the formula 


T _ f , VTTrf _ Vi? 4 + l* 
L J T ^2 R 


(4.9) 


For fixed i s , L decreases with decreasing N and as a result T decreases. When we include 
the first 1/iV correction the time of collapse is [5] 


yiTTg , _ rg r*(l + 3(1 + rg)) 

— r 4 6A^ 2 (1 + rg) 3 / 2 y/rQ — r 4 6A^ 2 (1 + r 4 )y/1 + — r A 

(4.10) 

By performing numerical integration of the above for several values of the parameter ro and 
some large but finite values of N, we see that the time of collapse is smaller for the 1/N 
corrected case. This means that, in the region of large N the time of collapse decreases as 
N decreases, with both the leading large N formula and the 1/N correction being consistent 
with this trend. However, as we saw above the time of collapse at n — 1 and n = 2 are 
larger than at n = oo. This means that there are one or more turning points in the time of 
collapse as a function of n. 

The deceleration effect that arises in the comparison of n — 1 and n — 2 with large 
n may have applications in cosmology. Deceleration mechanisms coming from DBI actions 
have been studied in the context of bulk causality in AdS/CFT :2I) E22J and applied in the 
problem of satisfying slow roll conditions in stringy inflation |23j . Here we see that the finite 
n effects result in a further deceleration in the region of small n. 
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We turn to the proper acceleration which is important in checking the validity of our 
action. Since the DBI action is valid when higher derivatives are small, it is natural to 
demand that the proper acceleration, should be small (see for example eh). The condition 
is 7 3 £ s d^R p hy S <1. In terms of the dimensionless variables it is 7 3 (<9 2 f) <C VN. If we want 
a trajectory with initial radius ro such that the proper acceleration always remains less than 
one through the collapse, then there is an upper bound on r 0 . This upper bound goes to 
infinity as N ^ 00 . We are already constrained by the condition that the spatial derivatives 
are small, to lie within the small or large r region, for finite-Ah For small ro we are in the 
matrix theory limit and things are well behaved. For large r 0 and r-large, the acceleration 
is under control, a ~ 1/r and the velocity will be close to zero. Interestingly, there will also 
be a valid large ro, relativistic regime. Consider for example the n — 1 case. The proper 
acceleration can be written as 

2r 3 -3 + 2s 2 + r 4 (s 2 - 2) 

a-- l + r 4 vT^2(r4(s 2 -4)-3) ' > 

For s ~ 1 and small r, this becomes 


2 r 3 

a ~- , 

3%/r^s 2 


(4.12) 


and y/1 — s 2 


can be found from the energy at the same limits, in which (El becomes 


Vl — s 2 




1 

M 175 ' 


(4.13) 


Therefore, we can identify a region where the proper acceleration is small by restricting it 
to be of order 1 /ro for example 


2 r 3 

a ~- 

3 


(2rg 




1 

?"o 


(4.14) 


This means that in regions where r ~ r 0 5 ^ 9 , we will have a relativistic limit described by the 
DBI, where stringy corrections can be neglected. This result also holds in the large -N limit. 
It will be interesting to develop a perturbative approximation which systematically includes 
stringy effects away from this region. 

Another quantity of interest is the effective mass squared E 2 — p 2 , where p = dC / ds is 
the radial conjugate momentum. It becomes negative for sufficiently large velocities. This 
includes the above regime of relativistic speeds and small radii. It is straightforward to see 
that if our collapsing configuration is considered as a source for spacetime gravity, this implies 
a violation of the dominant energy condition. In the context of the BFSS matrix model, it 
has been shown that for an action containing a background spacetime Gjj = rjjj + hu , in the 
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linearised approximation, linear couplings in the fluctuation hoi correspond to momentum in 
the X 1 direction [23|. The same argument can be developed here for the non-Abelian DBI. 
We couple a small fluctuation h 0r , which in classical geometry we can write as hoi = h 0r Xi 
for the unit sphere. We replace ay by cq/n. The action for DO-branes P3B is generalised 
from m by replacing R in A d t $i = A(i?)cq = -^cq with (R + ho r )- It is then clear that 
the variation with respect to R , which gives p, is the same as the variation with respect to 
hor, which gives T 0r . Hence, the dominant energy condition will be violated, since E < \p\ is 
equivalent to T 00 < T 0r . The violation of this condition by stringy D-brane matter can have 
profound consequences. For a discussion of possible consequences in cosmology see [25] . In 
this context, it is noteworthy that the violation can occur near a region of zero radius, which 
could be relevant to a near-big-bang region in a braneworld scenario. 


4.3 Distance to blow-up in D 1 _L D 3 


Comparisons between the finite and large N results can be made in the spatial case using 
the conserved pressure. The arguments are similar to what we used for the time of collapse 
using the energy functions. Consider the case n — 1, and let P = P/N. First calculate the 
derivative of the pressure - 


dP s(4r 4 + r 4 s 2 + 3) 
~dl ~ % /TT7 4 (l - s2)V2' 


This is clearly always positive. Now evaluate 


P 




(I + r o) 1/2 

1 + r 4 


(l + r 4 + r 4 ). 


(4.16) 


This should be compared with P(r, s ), evaluated for the value of s which solves the n — 1 
equation of motion, which by conservation of pressure is — ( 1+2r °) 


\A+hT 


Take the difference to find 


P 


r,s — 


\Jr A — ? 

y/iTVl) ' y/TTT 4 v / TT7 4 (i + r 4 )' 


+ 


(1 + 2r 4 ) 4(r 4 - rp) 


(4.17) 


Thus at fixed ro and r, P n= i, when evaluated for the value of s which solves the large 
n problem, is larger than when it is evaluated for the value of s which solves the n — oo 
problem. Since P increases monotonically with s for fixed r, this shows that for fixed r 0 , and 
any r, s is always larger in the large N problem. Since S = f dr/s, this means the distance 
to blow-up is smaller for n = oo. Hence for fixed ro, the distance to blow-up is larger at 
n — 1. 
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5 Towards a generalisation to higher even-dimensional fuzzy-spheres 


For generalisations to higher dimensional brane systems, and to higher dimensional fuzzy 
spheres 0 EH 1121 mu, it is of interest to derive an extension of the expressions for the 
symmetrised traces given above. In the general case, we define N(k,n ) to be the dimension 
of the irreducible representation of SO(2k+l) with Dynkin label ..., which contains 
k entries. We then take C(m, k,n) to be the action of the symmetrised trace on m pairs of 
matrices X l} where i — 1,..., 2k + 1 

C(m,k,n) = ——j — -STr [X i X i ) m . (5.1) 

A(/c,n) 

Finding an expression for C(m, k,n ) is non-trivial. Investigations based upon intuition from 
the ADHM construction lead us to conjecture that for n odd 


C(m, k, n ) 


2‘n‘ 1 .i( 2m ~ 1 + 2 »i) 


"+i 

zliKd+^-M) 2 )^- 1 ) 

*3 = 1 14=1 


(5.2) 


while for n even 1 


C(m, k, n ) 


2 fc n,i=i(2m 1 + 2 ii) |- y-j- i' I'n . 



(2 i„) 2m 


(5.3) 


We give the arguments leading to the expressions above in Appendix A. 

For higher even spheres there will be extra complications at finite-n. Consider the case of 
the fuzzy S' 4 for concreteness. The evaluation of the higher dimensional determinant in the 
corresponding non-Abclian brane action will give expressions with higher products of dt^i 
and dkj = [$*, Tj] 


S 


-T 0 JdtSTr |l + A 2 (d t ^) 2 + 2A 2 <I),/I) /; + 2A 4 (^^) 2 - 4A%$ ifc $ A + 

A 6 1 1/2 

+ 2A 4 (d t $i) 2 $ jk $ kj - 4A 4 d t $ l $ ij $ jk d t $ k + —(e ijk imdt$i$ j k$im) 2 > ■ (5.4) 


The ansatz for the transverse scalars will still be 


= R(t)Xi , 


where now i — 1,... ,5 and the A"*’s are given by the action of SO (5) gamma matrices on the 
totally symmetric ?r-fold tensor product of the basic spinor. After expanding the square root, 

1 For in = 0 the value STr{Xi Ah) 0 = 1 is once again imposed. 
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the symmetrisation procedure should take place over all the Xf s and [X l} Aj] 's. However, the 
commutators of commutators [[A", A"], [A, A]] will give a nontrivial contribution, as opposed 
to what happens in the large-n limit where they are sub-leading and are taken to be zero. 
Therefore, in order to uncover the full answer for the finite-n fuzzy S 4 it is not enough to just 
know the result of STr(X i X i ) m - we need to know the full STr ((A A) mi ([A, A][A, A ])) m2 
with all possible contractions among the above. It would be clearly interesting to have the 
full answer for the fuzzy S' 4 . A similar story will apply for the higher even-dimensional fuzzy 
spheres. 

Note, however, that for R = 0 in O all the commutator terms will vanish, since 
they scale like R 2 . This reduces the symmetrisation procedure to the one involving X t X t 
and yields only one sum for the energy. The same will hold for any even-dimensional S 2k , 
resulting in the following general expression 


£n,*(0,s) = -STrf^(-l) m S 2m (2m-l)(A i A i ) m [ 1/2 ') 

m =0 \m J 

= -#,n)^(—ir/ m (2m-l)C(m,fc,n)(^ 2 ) • 


m =0 


Using (Q, notice that in the odd n case 


2k—1 


n -\-1 
2 


C(m, k, n) = J] k, l)± 3 - i)*» 

*2 = 1 ' U + l2 ’ 


foddi 1, k, 1) 


The factor f nr i,i is 


k-1 


TT / / ^ \ 2 / 1\^ 

fodd{h,k,n) = [[11-+UJ - (U - - 


24 = 1 


(5.5) 


(5.6) 


(5.7) 


Inserting this form for C(m, k,n) in terms of C(m, k , 1) we see that 


2k—1 

£ n ,k{ 0, s) = N(n, k) JJ 
22 —1 


(1 + * 2 ) 

(n + i 2 ) 


n+1 

2 

E 

*3=1 


foddij' 3 1 ^5 
foddi 1 , k, 1 ) 


4=i,fc(0, s(2* 3 - 1)) 


(5.8) 


Similarly we derive, in the even n case, that 


where 


2 fc —1 


£„.»(0, s) = N(n, k) f[ AjA ! ° 


t (*3 ,k,n) 


• + A) • feven{ 1; k, 2) 

*2 — J- *3 — 




fc-1 

feven{.i 35^?^) 

24 = 1 




(5.9) 


(5.10) 
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and S is the energy density, i.e. the energy divided a factor of N(n, k). 
It is also possible to write C(m, k, 2) in terms of C(m, k, 1) 


C(m, k, 2) 


fc-i 

— 2 2m C(m, k, 1) n 

24=1 


24(24 + 2) 
^4(^4 +1) 


2 fc —1 


n 

22=1 


(^2 + 1) 
(*2 + 2) 


= 2 2m C(m, k, 1) 


fevenly 1 ; k, 2 ) 

/odd(l, k, 1) 


2fe—1 

n 

22 = 1 


(^2 + 1) 
(*2 + 2 ) 


which is valid for all values of m/0. 


(5.11) 


It turns out to be possible to give explicit forms for the energy for the n — 1 and n = 2 
case. Since the definition of the physical radius in section 2 is also valid for higher dimensional 
fuzzy spheres, we can express the results in terms of the rescaled variables r and s 

^71=1,fc(0) 

^n=2,fc(0;^) 

When plugged into (15.8|) . (15.91) the above results provide a closed form for the energy at 
f = 0, for any n and any k. 


2fc+l 

(1 - S-) 2 


(k + 1) 


(1 


2fc+l 


(2 k + 1) 


(5.12) 


6 Summary and Outlook 

We have given a detailed study of the finite N effects for the time dependent DO — D 2 fuzzy 
sphere system and the related D 1 _L D 3 funnel. This involved calculating symmetrised 
traces of SO(3) generators. The formulae have a surprising simplicity. 

The energy function E(r, s ) in the large N limit looks like a relativistic particle with 
position dependent mass. This relativistic nature is modified at finite N. Nevertheless our 
results are consistent with a fixed relativistic upper speed limit. This is guaranteed by an 
appropriate definition of the physical radius which relies on the properties of symmetrised 
traces of large numbers of generators. We showed that the exotic bounces found in the large 
N expansion in (Hj do not occur. It was previously clear that these exotic bounces happened 
near the regime where the 1/iV expansion was breaking down. The presence or absence of 
these could only be settled by a finite N treatment, which we have provided in this paper. 
We also compared the time of collapse of the finite N system with that of the large N system 
and found a finite N deceleration effect for the first small values of N. The modified E(r, s ) 
relation allows us to define an effective squared mass which depends on both r, s. For certain 
regions in (r, s) space, it can be negative. When the DO — D 2 system is viewed as a source 
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for gravity, a negative sign of this effective mass squared indicates that the brane acts as a 
gravitational source which violates the dominant energy condition. 

We extended some of our discussion to the case of higher even fuzzy spheres with SO(2k + 
1) symmetry. The results for symmetrised traces that we obtain can be used in a proposed 
calculation of charges in the D1 _!_ D(2k + 1) system. They also provide further illustrations 
of how the correct definition of physical radius using symmetrised traces of large powers of 
Lie algebra generators gives consistency with a constant speed of light. A more complete 
discussion of the finite N effects for the higher fuzzy spheres could start from these results. 
Generalisations of the finite N considerations to fuzzy spheres in more general backgrounds 
[2> will be interesting to consider, with a view to possible applications in cosmology. 
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Ward for useful discussions. The work of SR is supported by a PPARC Advanced Fellowship. 
SM and CP would like to acknowledge a QMUL Research Studentship. This work was in 
part supported by the EC Marie Curie Research Training Network MRTN-CT-2004-512194 
Superstrings , and by the PPARC Research Grant PPA/G/0/2002/00477 M Theory, String 
Theory and Duality. 
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A General formula for the symmetrised trace 


As in section 5, we define N(k,n ) to be the dimension of the irreducible representation of 
SO(2k + l) with Dynkin label (|,|,...,|). These are the usual fuzzy sphere representations 
0 m (f° r example, for k=l the X t are the elements of the Lie algebra of SU(2) in the 
irreducible representation with spin (|). Then 


N(k,n) 


n 


n + 2k — (i + j) + 1 
2k - (i — j) — 1 


-A- n + 2k ~ 21 + 1 

2k - 2/ - 1 
z=i 


(A.l) 


The symmetrised trace is dehned to be the normalised sum over permutations of the 
matrices 

STr(X h ... X ip ) = i y Tr(X, m ... X„ w ) . (A.2) 

P ' aeSp 


We have given earlier a conjecture for the symmetrised trace of m powers of the quadratic 
Casimir A'A* = c I nXre , where c = n(n + 2k). This is, for all m, k and n even 


1 

N(k,n) 


STr(X i X l ) m 


2 *' nf;_i(2m -1 + 2 »q 


f k-1 

*3 = 1 *4 = 1 



(A.3) 


where for k — 1 the product over i 4 = 1,..., k — 1 is just dehned to be equal to 1. Similarly 
for all m, k and for n odd we have proposed that 


1 

N(k,n) 


STr{X i X l ) m 


2 k Yl- 1= i(2m - 1 + 2t,) 

(fc-i)in ?=iV+* 2 ) 


£[n((y-) 

*3 = 1 *4 = 1 


(i 3 -i) 2 )(2! 3 -l) 2m 


(A.4) 

The argument leading to these conjectures follows. Firstly, it is possible to view the fuzzy 
sphere matrices X t as the transverse coordinates of the world-volume theory of a stack of D 1- 
branes expanding into a stack of D(2k + l)-branes [3 27, T8j. There is also a dual realisation 
of this system in which the Dl-branes appear as a monopole in the world-volume theory of 
the D(2k + l)-branes. The ADHM construction can be used to construct the monopole dual 
to the fuzzy sphere transverse coordinates. If one takes the N(k, n)-dimensional fuzzy sphere 
matrices representing a stack of N[k , n) Dl-branes as ADHM data for a monopole, then one 
naturally constructs a monopole dehned on a stack of N{k — 1, n + 1) D(2k + l)-branes. We 
have also done the calculation which shows that the charge for the monopole just constructed 
gives precisely N(k,n ), which provides a consistency check (this calculation extends some 
results in and will appear in EDI). 


It is also possible to calculate the number of D(2k+1) branes from the fuzzy sphere ansatz 
for the transverse coordinates, by looking at the RR coupling on the D-string worldvolume. 
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In this case one does not get N(k — 1, n + 1) as one would expect, but instead the quantity 

^_ W l)Il£lV±i), (A.5) 

C 2 

where c = n{n + 2k). This number of branes does agree with N(k — 1, n + 1) for the first 
two orders in the large n expansion 

Number of D(2k + l)-branes = N{k — 1, n + 1) ^1 + O • (A.6) 


Now consider this RR charge calculation more carefully. First take the k= 1 case. Based 
on the ADHM construction we expect the number of H3-branes to be lV(0,n + 1) = 1. 
However, equation (H3 suggests that the charge calculation gives for k = 1 the answer 

n + 1 

T 

C2 

Suppose that the numerator in the above is correct, but that the denominator is correct only 
at large n and that it receives corrections at lower order to make the number of Z)3-branes 
exactly one. Then these corrections need to satisfy 


(A.7) 


13 5 

1 = (n + l)(c“2 + rzqcT 2 -|- x 2 c~ 5 + ....). (A.8) 

It is easy to show that we need x\ — — \ and x 2 = §, by Taylor expanding and using that 
c = n{n + 2) for k — 1. Therefore, we would like to have a group theoretic justification for 
the series 


_i 1 3 5 

c 2 - -c 2 + -c 2 + ... . 

2 o 

There exists a formula for the first three terms in the large n expansion of the k 
symmetrised trace operator pi], namely 


(A.9) 
= 1 


—- -STr(XiXi) m = c m — m(m — l)c m H- m(m — l)(m— 2)(7m — l)c m +... . (A.10) 

N(l,n) 3 45 


Now, if we make the choice m = —\ in (1A. 101) we get precisely (IA.9I) . However, this suggests 
that if this choice is correct, then we should have an all orders prediction for the action of 
the symmetrised trace operator. Thus, for k = 1 we predict that 


N(l,n) 


STr(X i X i ) r 


(n + 1) ’ 


(A.11) 


where for future reference we consider the left hand side to be equal to the symmetrised 
trace in a large-n series expansion, as appeared in [5]. 

Checking the conjecture (1A. 1 If) beyond the first three terms in a straightforward fashion, 
by techniques similar to those employed in |5|, proves difficult. This involves either adding 
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up a large number of chord diagrams, or complicated combinatorics if one uses the highest 
weight method. 


An alternative approach involves first writing down the conjecture based on brane count¬ 


ing for general k, since the methods of [5] turn out to generalise from the k — 1 to the general 
k case. The conjecture for general k, based on the brane counting, follow immediately from 




(A. 12) 


Note that the right hand side of this equation appears in the factor outside the sum in (D 
and (IA.4j) . Notice also that the above expression concerns the large n expansion of the 
symmetrised trace considered at m — — k + 


One can repeat the k — 1 calculation of p] for general k, to check the first three terms 
of this conjecture. A sketch of this calculation follows before displaying the full results. 


First we calculate < STr{XiXi) m for m = 2,3,4. Then we find the first three terms in 
the symmetrised trace, large n expansion using these results. Finally we can check that 


the conjecture ET1 is true for the first three terms in the symmetrised trace large n 


expansion, for general k as well as for k — 1. We then proceed to calculate the fourth 
term in the expansion, for general k. To do this we need to calculate n - ) STr(XiXi) m for 

m = 5,6. We then show that the fourth term in the large n expansion of the symmetrised 
trace agrees precisely with (IA.12I) . 

In the following, we use the notation of j5], with each trace of a string of 2m A 7 , matrices 
arising here being represented by a chord diagram with m chords. This provides a convenient 
way to represent equivalent strings of matrices. 

For the calculation of STr(AjA*) 2 there are three different strings of the four X t matrices 
and two different chord diagrams. Two of the three strings correspond to the same chord 
diagram. In the following, the first column contains a fraction which is the multiplicity of 
the chord diagram in the list of strings divided by the total number of strings. The second 
column contains a picture of the chord diagram preceded by an example of a string in the 
equivalence class defined by this chord diagram. The evaluation of the chord diagram is the 
final entry. 


I 1212 


| 1122 
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Using this, one finds immediately that 


—i-^S7Y(A'.A'i) 2 = c 2 - \kc . 
N(k,n) 3 


(A.13) 


For m — 3 there are 15 different strings of X t matrices and five different chord diagrams, 
which evaluate as follows 



Thus we find that 

1 16 

— STr^Xif = c 3 - 4 kc 2 + — k{Ak + l)c . (A.14) 

J\{k:, n) 15 


For rn — 4 there are 105 different strings and 18 different chord diagrams. We omit the details 
for simplicity. The final result is that 


N(k, n 


I 1^ 

-STrfXiXj) 4 = c 4 - 8 kc 3 + —k(7k + 2)c 2 - 

5 


105 


k(34k 2 + 24k + 5)c. (A.15) 


For rn — 5 there are 945 different strings of Xj matrices, and 105 different chord diagrams 2 . 
The result is 


1 40 16 64 

—- -STriXiXif = c 5 - kc 4 +—(13k + A)kc 3 - (158k 2 + 126k + 31)kc 2 

N(k,n) 3 3 63 


256 

+ — (496U 3 + 672k 2 + 344A; + 63)A;c. 


(A.16) 


2 We acknowledge the assistance of Simon Nickerson, for writing a computer programme used here. Maple 


files for these calculations are available from the authors. 
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For m = 6 there are 10395 different strings of X % matrices, and 902 different chord diagrams, 
and we find that 


16 4 64 2 


N(k, n 


-STr(XiXi) b = c b — 20/cc 5 H- (31k + 10)fcc 4 — —(677k 2 + 582A; + 157)kc' 


256 


315 
1024 

10395 


63 


(1726fc 3 + 2616fc 2 + 1541fc + 336)fcc 2 


(11056/c 4 + 24256/c 3 + 22046& 2 + 9476A: + 1575)fcc . 


(A.17) 


Now we calculate the first four terms in the large n expansion of STr(XiXi) m . Suppose 
that the coefficient of c m ~ l term in STr(XiXi) m is a polynomial in m of order 21. Then we 
have the following ansatz: The known factors of these polynomials come from the fact that 
the series has to terminate so that there are never negative powers of c for m = 1,2,3,... 
Then 


N(k,n 


■STr(X i X i ) m = c m + y 4 (k)m(m — 1)6 


rn— 1 


m + V 3 {k) \m(m — 1 )(m — 2)6 


,771—2 


m 2 + y 5 (k)m + y§(k) )m(m — 1 )(m — 2)(m — 3)c : 


m— 3 


+0(6 


, 771—4 \ 


We now find the unknown functions yi(k),y 2 (k).. .y§(k) using the results of STr(XiXi) m 
for m = 2, 3,4, 5, 6 calculated above. We find that 

Vi(k) = ~k, y 2 (k) = ^(5k + 2)k, 

m(k) = ^(k-2)k, y 4 (k) = ^^(-140/c 2 - 168fc - 64)fc, 

y»(*0 = ^(-84t 2 + 216t+192)i, y e (k) = ^-F(128fc 2 + 96k - 104)fc . 

(A.18) 


Now we are able to provide a check of the conjecture 1A .1 21) . First we express the right- 
hand side of (IA.12I) as a function of c as 


2k—1 

n 


L = -^ n + l y/(k 2 + c) A1 c + 2kl- l 2 


k -1 

n 


= C~ k+ 2 


oo 7 
c j ’ 


3 =0 


(A.19) 
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where 



b 2 - -(Wk 2 -Zk + 2)k(k--)(k + -)(k + ^) : 

63 = ——(—24 + 34k — 61k 2 + 56k 3 — 140fc 4 ) x 
2835 

k ( k ~ 5 ) i k + 5 ) i k + 1 ) i k + 5 ) ■ ( a - 2o) 


Now consider the left-hand side of (1A. 121) involving the large n expansion of STr(A 7 jA 7 j) m , 
which we calculated above, but now we set m — — k + Expanding in inverse powers of c, 
we find that this becomes 


N(k, n 


-STriXiXiY 


m=-k-\- 1/2 


= C~ k+ 2 


3 V- 

c j 


3=0 


with precisely the coefficients b{ given in (IA.20I) . Given the extensive and non-trivial calcu¬ 
lations required to obtain these results, we believe that there is strong evidence for the truth 

of ira . 

For k — 1 the guess of the exact answer for n even is 


—L-STr(X t X,r = 2(2 ”+ 1) Z< 2 if m ■ (A.21) 

It is easy to show that (IA.21I) agrees with the first four orders in the large n expansion (JA.1 81) 
for k — 1. If we set m = —| in this we get zero, because of the (2m + 1) factor. This might 
appear to contradict (1A.11I) . but it is easy to show, using a large n expansion, that if (1A. 211) 
is true then (IA.11D holds to all orders. To calculate the large n expansion of this sum we can 
use the Euler-Maclaurin formula. This approximates the sum by an integral, plus an infinite 
series of corrections involving the Bernoulli numbers i? 2 P 


E/w 


fn+l 


f(x)dx + ~[f(n+ 1) — /(O)] 

00 JD 

+ E(^fT[/ (2P " 1) ('*+ 1 )-/ (2p_1) (0)]' 


(A.22) 


We see from this calculation that for k — 1 the value m = — | is very special. It is the only 
value of m for which the higher order terms in the Euler-MacLaurin large n approximation 
of the sum in (IA.21I) are zero. 
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B Finite n results on symmetrised traces from the highest weight 
method 


Results on finite n symmetrised traces can be obtained by generalising the highest weight 
method of j)|. For the 50(3) representations used in fuzzy 2-spheres we have 

i5Tr J= 1 / 2 (aiai) m = (2m + 1) , (B.l) 

where the 1/2 comes from dividing with the dimension of the spin- 1/2 representation. A 
similar factor will appear in all of the results below. The above result was derived in jo). For 
the spin one case, we will obtain 


lorn , srn 2 2m+1 (2m + 1) 

-STr J=1 («;«;) =--- 


(B.2) 


These results can be generalised to representations of 50(2/ + 1) relevant for higher fuzzy 
spheres. The construction of higher dimensional fuzzy spheres uses irreducible representa¬ 
tions of highest weight (^,••• , |), as we have noted. For the minimal representation with 
n = lwe have 

(2/ + 2m — 1)!! 


-5Tr n= 1 (WW) = 


(B.3) 


D n =i -'." (2m — 1)!!(2/ — 1)!! 

Notice the interesting symmetry under the exchange of l and m. For the next-to-minimal 
irreducible representation with n — 2 we obtain: 


1 

D n =2 


STr n=2 (XjXi) 


2 2m (l + 1) 


(2/ + 2m — 1)!! 
(2m- 1)!!(2/ + 1)!! ' 


(B.4) 


This is a generalisation of the spin 1 case to higher orthogonal groups. It agrees with the 
formulae in section A of the appendix, with / —> k. 


B.l Review of spin half for 50(3) 


We will begin by recalling some facts about the derivation of the n 
commutation relations can be expressed in terms of 03 , a± 


a± 


[o 3 , a±] 


o_] 

c 


1 

7! 


(«i ± ia 2 ), 


2a±, 

2o 3 , 


T cc_q: + T ct 2 . 


1 case in j3]. The 


(B.5) 
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With these normalisations, the eigenvalues of a 3 in the spin half representation are ±1 and 
a + a_ is 1 on the highest weight state. 

It is useful to define a quantity C(p, q) which depends on two natural numbers p,q and 
counts the number of ways of separating p identical objects into q parts 


C(p, q) 


(P + g-1)! 

p\(q- 1)! 


(B.6) 


We begin by a review of the spin half case, establishing a counting which will be used again 
in more complicated cases below. This relies on a sum 


2n—2k is 12 

2t E -"EE*- 1 ) 

hk * 2 = 011=0 


*l+*2+---*2fc 




Recall that this sum was obtained by evaluating a sequence of generators of SO (3) consisting 
of k pairs a_a + and with powers of a 3 between these pairs - 


a\ J2k+1 a + a\ Jk a- ■ ■ ■ a-a J z 3 a+a^ 2 a-a^ 1 . 


(B.7) 


We can move the powers of a 3 to the left to get factors (a 3 — 2) j2+j4+ '" j2fc . Moving the a 3 
with powers Ji, J 3 .. gives a^ 1+Ji+ '". The k powers of a_a + gives 2 k . The above sum can be 
rewritten 


2m—2k 2m—2k—Js~\--.J2k+i 2m— 2k— J 2 ... J 2 k+i 


E 

+fc+l=0 


E 

e/ 2 — 0 


E 

Ji=0 


- 1 ) 


Ji+Ji-'-Jik — 2 k 


n\ 


(n — k)\k\ 


(B.8) 


This includes a sum over J e = J 2 + J\ + ... + J 2 k■ The summand does not depend on 
the individual J 2 , J 4 ,... only on the sum J e which ranges from 0 to 2 m — 2k. The sum 
over J 2 , J 4 ,... is the combinatoric factor, introduced above, which is the number of ways of 
splitting J e identical objects into k parts, i.e. C(J e , k). The remaining 2m — 2k — J e powers 
of a 3 are distributed in k + 1 slots in C(2m — 2k — J e , k + 1) ways. Hence the sum (IB. 811 can 
be written more simply as 


2m— 2k j 

2 k V (-1 ) J *C(J ei k)C(2m - 2k - J e , k + 1) = 2 k - ——ttttt • (B.9) 

t—* (m — ky.kl 

J e —0 


Then there is a sum over k from 0 to m, with weight 


C(k,m) 


2 k k\(2m — 2 k)\m\ 
(■m — k)\(2m)\ 


(B.10) 


which gives the final result 2m + 10. Similar sums arise in the proofs below. In some cases, 
closed formulae for the sums are obtained experimentally. 
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B.2 Derivation of symmetrised trace for minimal SO(21 + 1) representation 


The Casimir of interest here is 

X,x, = Xi l+1 + ^ . (B.ll) 

i=1 


The patterns are similar to those above, with 03 replaced by X 21 + 1 , and noting that here 
there are l “colours” of a± which are X± . All the states in the fundamental spinor are 
obtained by acting on a vacuum which is annihilated by l species of fermions. Generally we 
might expect patterns 

V * 1 vCh) Vl2 vlh) /T 3 1 Q\ 

•■■ A 2Z+1 A - ^2Z+1^+ ••• (D.lZ) 

In evaluating these, we can commute all the X 2 j+i to the left. This results in shifts which 
do not depend on the value of j. It is easy to see that whenever X^ is followed by X^ 
we get zero because of the fermionic construction of the gamma matrices. cannot also 
be followed by X+ 1 because X+ - 1 X+ ! + X+ - 1 X+ 1 = 0. So the pairs have to take the form 
X^X^ for fixed j. The sum we have to evaluate is 


m 2m—2k 

EE (-1 ) Js C(J e , k)C(2m -2 k- J e , k)C(k, l)2 k C(k, m ) 

k =0 J e = 0 

m , \ 

= J 2‘C(t,m)C(k,l) (B.13) 

k =0 ^ ' 

(21 + 2m - 1 )!! 

“ (2m- 1)!!(2Z - 1)!! ' 

The factors C(J e , k) and C(2m — 2k — J e , k) have the same origin as in the spin half case. 
The factor (IB. 101) is now generalised to an Z-colour version 

C(ki,k 2 ...ki;m) = 2 k — — m '^ ki\k 2 \...ki\ (B.14) 

(2m )\ (m — k)\ 

This has to be summed over k\,ki. For fixed k = k\ + ■ ■ ■ + ki we have 

u 

V C(k 1 ..k h m)——— = V C(k, m) = C(k , m)C(k, l ) . (B.15) 

' ky.-ky. 

kt-ki ki..ki 

The combinatoric factor , in the second line above conies from the different ways of 
distributing the k\..kj pairs of (—h) operators in the k positions along the line of operators. 
The subsequent sum amounts to calculating the number of ways of separating k objects into 
l parts which is given by C(k, l). The C(k, m ) is familiar from (IB.101) . This sum can be done 
for various values of k, m and gives agreement with (IB.31) . 
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B.3 Derivation of spin one symmetrised trace for SO(3) 

For the spin one case more patterns will arise. After an cm acts on the highest weight, we 
get a state with 03 = 0 so that we have, for any positive r 

a^a-\J = 1,03 = 2 >= 0 , V r > 0 . (B.16) 

Hence any cm can be followed immediately by a + . These neutral pairs of (aq_cm) can be 
separated by powers of 03 . Alternatively an cm can be followed immediately by cm. The 
effect of a 2 _ is to change the highest weight state to a lowest weight state. I 11 describing the 
patterns we have written the “vacuum changing operator” on the second line, with the first 
line containing only neutral pairs separated by 03 ’s. Let there be J\ neutral pairs in this 
first line and L\ powers of a 3 distributed between them. After the change of vacuum, we 
can have a sequence of (cmaq_) separated by powers of a 3 . Let there be a total of J 2 neutral 
pairs and L 2 c^’s in the second line. At the beginning of the third line we have another 
vacuum changing operator a 2 + which takes us back to the highest weight state. In the third 
line, we have J 3 neutral pairs and L 3 powers of a 3 . The equation below describes a general 
pattern with p pairs of vacuum changing operators. The total number of neutral pairs is 
2 p + J where J = J\ + J 2 + • • • + J 2p+ ±. The general pattern of operators acting on the 
vacuum is 

# (a+cm) # ( a + a _) # ■ ■ ■ # (a+cm) # \J — 1, a 3 — 2 > 

# (cm«+) # (cma+) # ■ ■ ■ # (cma + ) #cr 
# («+«_) # ( a + a _) # • • • # (a+cm) # a 2 + 

# (cma+) # ( a-a + ) # • • • # (cma + ) #a 2 _ 

# (a+a-) # (a+cm) #■ • ■# (a + a_)# , (B.17) 

where in the above the first line of operators acts on the state \J = l,a 3 = 2 > first, 
then the second line acts, and so on. The symbols # represent powers of a 3 . We define 
J e = J 2 + J 4 + • • • J 2p which is the total number of (—h) pairs on the even lines above. 
There is a combinatoric factor C(J e ,p) for distributing J e among the p entries, and a similar 
C(J — J e ,p + 1) for the odd lines. The L e = L 2 + L 4 + • • • + L 2p copies of a 3 can sit in 
(J 2 + 1) + (J 4 + 1) + • • • + (J 2p + 1) positions which gives a factor of C(L e , J e + p). The 
L\ + L 3 + • • • + L 2 p+i can sit in (Ji +1) + ( J 3 +1) + • • • + (J 2p +i + 1) = J — J e -\-p-\- \ positions, 
giving a factor C(2m — 2J — Ap — L e , J — J e + p + 1). There is finally a factor C(2p + J, m) 
defined in (IB. 1 0il which arises from the number of different ways the permutations of 2m 
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indices can be specialised to yield a fixed pattern of a + , cm, «3 


[m/2] m—2p J 2m—Ap—2J 

EEE E c( j e , P ) c(j-j e , P +i)(-i) L ' c(L',j e + P )x 

p=0 J= 0 J e = 0 L e =0 

C(2m — 2 J — Ap — L ei J — J e + p + 1) x 

2 2m-2J-4 P Q( 2j i)J« Q(2, 2) p C(2p + J, m) . 

By doing the sums (using Maple for example) for various values of m we find -- &m+i) . 

The factors denoted in |S] by N(i,j), arise from evaluating the cm, a+ on the highest 

weight. 

B.4 Derivation of next-to-minimal representation for SO(21 + 1) 


The n = 2, general l patterns are again similar to the n — 2, l — 1 case except that the 
o_,o + are replaced by coloured objects of l colours, i.e. the X±\ We also have the simple 
replacement of a 3 by X 2 i+\. 

We define linear combinations of the gamma matrices which are simply related to a set 
of / fermionic oscillators: T 1 / = ^(T 2 j_i +iT 2 j) = \/2a\ and T^ = ^(T 2 j_! — iV 2 i) = \[2ai. 
As usual X t are expressed as operators acting on an n-fold tensor product, and 



(B.18) 


Some useful facts are 


x r 2 l+ 1 x + 

0 

> 

X 2l+1 X_X + 

0 

> 

x_x + x r 2l+1 xl 

0 

> 

X+Y+X+ 

0 

> 

X+X_Xl 

0 

> 

x'ix'l 

0 

> 


0, X r 2l+1 Xl\0> = (-2) r Xl\0>, 

X_X + X r 2l+l \0 > = (2) r X_X + \0 >, 

0, Y + Xl + XlY + \0> = 0, 

0, X.Y+X+ 10> = 0, 

Q( 2,l)Xj|0>, X+X_X+Y+\0 > = Q(2,l)X+y+|0>, 

Q(2, 2)10 >, Y_X_X + Y + 10> = Q(2,2)|0> . 


It is significant that the same Q(2, 1), Q(2, 2) factors appear in the different places in the 
above equation. In the above X + stands for any of the l x/’s. Any equation containing 
X± and Y± stands for any pair X± and X^P for i,j distinct integers from 1 to l. 

The general pattern is similar to (IB. 171) with the only difference that the (cc_a + ) on the 
first line is replaced by any one (X_ x/ ) for i = 1,...,/. The positive vacuum changing 






operators can be (X^X^), where i, j can be identical or different. For every such choice the 
allowed neutral pairs following them are X+ X^ and the dual vacuum changing operator is 

The summation we have to do is: 


[m/2] m—2p 2m—4p— 2 J J , 

EE E E C(2p + J, m)C(2p + J, l)C(J e , p) C(J -J e ,p+ 1) x 

p=0 J= 0 L e = 0 Je= 0 

(— l) Le C(L e , J e + p)C(2m — 2 J — Ap — L e , J — J e + p + 1) x 
2 2m—2J— 4p Q^iy-Je g( 2 ,l) Je g( 2 , 2 ) p ^ . 

The g-factors can be easily evaluated on the highest weight and then inserted into the 
above 

g(l,l)=4, g(2,l) = 4 g(2,2) = 16. (B.19) 

By computing this for several values of m, l, we obtain (El- Note that both the l = 1 and 
the general l case will yield the correct value for m — 0, which is 1 . 
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